Subdiffusive axial transport of granular materials in a long drum mixer 
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Granular mixtures rapidly segregate radially by size when tumbled in a partially filled horizontal 
drum. The smaller component moves toward the axis of rotation and forms a buried core, which 
then splits into axial bands. Models have generally assumed that the axial segregation is opposed by 
diffusion. Using narrow pulses of the smaller component as initial conditions, we have characterized 
axial transport in the core. We find that the axial advance of the segregated core is well described by 
a self-similar concentration profile whose width scales as t^ , with a ^ 0.3 < 1/2. Thus, the process 
is subdiffusive rather than diffusive as previously assumed. We find that a is nearly independent 
of the grain type and drum rotation rate within the smoothly streaming regime. We compare our 
results to two one-dimensional PDE models which contain self-similarity and subdiffusion; a linear 
fractional diffusion model and the nonlinear porous medium equation. 

PACS numbers: 46. lO.+z, 64.75. +g 



An interesting property of dry granular materials is 
their tendency to s^arate by size under a wide variety 
of flow conditions 1,*^. Granular segregation is widely 
found in nature, and plagues industrial processes as well. 
Probably the best controlled and most widely studied ex- 
ample is segregation along the axis of a partiall y fi lled, 
horizontal "drum mixer" 3, 4, 5, |1 0, |i |1 EHH El 
E^Q^Q- After hundreds of drum rotations, an initially 
mixed binary distribution of different-sized grains sorts 
itself into almost periodic bands along the axis of the 
drum. These bands are threaded by a radial core of the 
smaller grains which develops prior to axial band forma- 
tion J,, J^, JjJ, Q, [141 • The radial core typically forms 
after just a few drum rotations. Accounting for this rich 
dynamical behaviour has been the goal of cellular au- 
tomata models 16], molecular-dynamics simulations 1 181 
and several continuum theories llllllllil,!!!!!!!!!^^^ 
The axial bands must somehow be sustained against be- 
ing mixed away by the random motion of the grains. 
Continuum models have generally assumed that the ran- 
dom motions mimic normal diffusion, and therefore that 
normal Laplacian gradient terms determine the short- 
wavelength cutoff of the axial band pattern. In this Let- 
ter, we experimentally challenge this common assump- 
tion. We find, surprisingly, that the axial transport of the 
radially segregated core along the drum is much slower 
than diffusion, i.e. that it is subdiffusive. It is nev- 
ertheless described by a self-similar profile which scales 
approximately as t^^^. We also find that the self-diffusion 
of the larger particles is subdiffusive. These results have 
strong implications for models of axial segregation, and 
possibly for other theories of granular mixing. 

Early theoretical models regarded axial band forma- 
tion as the result of a diffusion process with a negative 
diffusion coefficient These models ignore the 

radially segregated core, and they cannot account for the 
oscillatory transient that precedes axial band formation 
in some mixtures 0, 0, E|- This oscillatory travelling 



wave state apparently demands that the basic dynamics 
be at least second order in time. A later model due to 
Aranson et. al 23, reproduced both axial segregation 
as well as the oscillatory transient, while still ignoring 
the core. We have recently shown experimentally that 
this model is also inadequate 10], and that the core dy- 
namics are themselves oscillatory. Another model due to 
Elperin et. al [2^ regards axial segregation as resulting 
from a radial core instability leading to a spatially peri- 
odic thickening of the core. This model, unfortunately, 
also cannot account for the travelling wave state. In all 
cases, these models explain the short wavelength cutoff of 
the axial band pattern as the result of the supposed axial 
diffusion of the smaller grains. We show experimentally 
below that axial transport is not well-described by nor- 
mal diffusion. This is true of either the smaller grains in 
a binary mixture or of the larger grains in a self-mixing 
process. This falsifies a common, basic assumption of 
segregation models. 

A few studies have inve stig ated the axial transport of 
grains experimentally jl3, llJ, ll^ ll3 , but none have sys- 
tematically investigated the effects of varying grain type 
and drum rotation rate. Here we report experiments 
which characterize the axial transport of radially segre- 
gated grains using several different grain types and drum 
rotation rates, starting with a narrow pulse initial condi- 
tion. 

The drum mixer used in all experiments consisted of 
a horizontal Pyrex tube, 600 mm long with an inner di- 
ameter of 28.5 mm, rotated about its long axis at a con- 
stant rotation rate of 0.31 rev/s or 0.62 rev/s. The larger 
grains were either cubic white table salt or transparent 
glass spheres and had a size range of 300-420 /im. The 
smaller grains were either irregularly shaped black hobby 
sand or bronze spheres, with a size range of 177-212 /im. 
The filled volume fraction of the drum was 28 %. In or- 
der to reproducably fill the drum, the grains were loaded 
into a long U-shaped channel, which was inserted length- 
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wise into the drum and rotated to deposit its contents. 
This procedure ensures a uniform fihing fraction of the 
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FIG. 1: (Color online) a) An image of the shadow of the 
radial core formed by 177-212 /xm sand grains within 300-420 
/im salt grains, b) The detected edge used to determine the 
vertical extent h{x,t) of the core, c) The x integral of h is not 
constant in time, and is thus not proportional to the volume of 
small grains contained within the radial core, d) The integral 
of is constant in time and is proportional to the desired 
concentration. 



drum. To obtain reproducible, quantitative dynamical 
information, we used a pulse initial condition. The pulse 
was made by placing thin spacers in the U-shaped chan- 
nel 1.5 mm apart. The 1.5 mm space was filled with the 
smaller grains, and the remaining space was filled with 
the larger grains. 

After a few drum revolutions, the pulse of small grains 
forms a subsurface radial core and cannot be observed 
using standard surface lighting and video imaging tech- 
niques. Instead, we used a bulk visualisation technique 
developed by Khan et al 0|. The large grains are 
transluscent and the small grains are opaque. When a 
bright light source is placed behind the rotating drum, 
one can observe a shadow of the radial core on the front 
face of the granular sample. This shadow is a two- 
dimensional projection of the radial core. A computer 
controlled high speed camera was used to observe the ra- 
dial core shadow. Five images per drum revolution were 
obtained and averaged to determine the evolution of the 
radial core. Figure^ shows a typical image of the radial 
core shadow. Using edge detection, the radial core 
height h{x^t) was measured as shown in figure Qd, and 
expressed as a fraction of the full height of the material 
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FIG. 2: (Color online) a) Concentration profiles of a spreading 
radial core pulse of sand grains within salt grains at various 
times, b) Power-law scaling of the FWHM of the radial core 
pulse. From the linear fit (yellow line) we find that the width 
oc t°"^®. c) Collapsed concentration profiles of the radial core 
pulse corresponding to (a). The collapse parameter is a = 
0.37. 



TABLE I: Collapse parameters for the self-similar spreading 
of radial cores in various grain types and rotation frequencies. 



Large grains 
300-420 /im 


Small grains 
177-212 /xm 


Rotation rate 
(rev/s) 
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salt 


sand 


0.31 


0.38 ± 0.03 


salt 


sand 


0.62 


0.37 =b 0.03 


glass 


bronze 


0.31 


0.31 =b 0.04 


glass 


bronze 


0.62 


0.29 =b 0.01 


glass 


sand 


0.31 


0.35 =b 0.03 


glass 




0.31 


0.34 ± 0.04 


salt 




0.31 


0.29 =b 0.01 



in the drum. If we assume that any cross section of the 
three dimensional structure of the radial core is an ellip- 
soid, the square of the diameter of the radial core h'^ is 
proportional to the volume of small grains contained in 
the radial core at each axial position x. Figure shows 



the time evolution of the x integral of /i, which increases 
with time. Figure shows that the x integral of h'^ is 
constant in time, as it should be for a conserved quantity. 
This validates our assumption about the core shape, and 
demonstrates that can be used as a local concentration 



measure. The error in measurement of h'^ corresponding 
to an error in of ± 2 pixels. 

For a normal diffusive process, the width of a nar- 
row pulse initial condition grows as t^/^. In our experi- 
ment, the pulse of small grains do not mix into the larger 
ones, but instead the pulse sinks below the surface of the 
larger grains forming a radial core, which then spreads 
axially. We can nevertheless ask if this axial spread- 
ing is analogous to normal diffusion, as is assumed in 
models [H Ea HI HI IH . Figure shows the radial 
core concentration profile at different times, for a mixture 
of small sand grains and large salt grains. Plotting the 
full-width at half-maximum of the concentration profile 
against time, we determined the power-law dependence 
of the radial core width with time, as shown in figure 
|21d. From this, we determined that the width scales as 
t", where a < 1/2. This analysis, however, only deter- 
mines the power-law time dependence of one arbitrar- 
ily chosen dimension of a pulse (here, the half-maximum 
width) and not the whole pulse shape. For a symmetric 
initial condition, data collapse can test the scaling of the 
entire pulse. Figure Et shows collapsed data correspond- 
ing to the concentration profiles in Et, where the axial 
length scale was transformed as X xt and the ax- 
ial concentration of small grains C{x,t) was transformed 
as C ^ Ct^. The pulse width increases at the same 
rate as the pulse amplitude decreases, thus the spreading 
process is self-similar. This implies that the integrated 
concentration is constant and that no grains are lost from 
the core. The average collapse parameter for large salt 
grains and small sand grains with a drum rotation rate 
of 0.62 rev/s is a = 0.37 ± 0.03, averaged over 10 runs. 
Similar experiments were repeated for different combina- 
tions of grains at two drum rotation rates. The results 
are shown in table 1 ffl We conclude that cores of small 
grains spread axially as where a ~ 1/3 < 1/2, inde- 
pendent of grain type and drum rotation rate within the 
smoothly streaming regime. 

It is interesting to compare the spreading of radially 
segregated cores of small grains with the non-segregating 
self-diffusion of the large grains alone. To observe this 
experimentally, some of the large grains were dyed black. 
These dyed grains were loaded into a drum full of other- 
wise identical white grains with a 1.5 mm wide pulse as 
the initial condition. The space-time evolution was ob- 
served using standard surface-lighting and imaging tech- 
niques H, 0- Figure int shows the concentration profile 
of dyed salt particles at various times. This data was col- 
lapsed in a similar way as discussed previously, as shown 
in figureint). Again, we find a collapse parameter a < 1/2. 
For runs using salt grains, a = 29 ±0.01 and for runs us- 
ing glass spheres, a = 34 ± 0.04, each averaged over 5 
runs. Thus, we conclude that the self-diffusion of grains 
in the rotating drum is also subdiffusive, even when no 
segregation is involved. We discuss some differences be- 
tween these two cases below. 




FIG. 3: (Color online) a) Concentration profiles of a mixing 
pulse of dyed black salt grains surrounded by white salt grains, 
b) Collapsed concentration profiles corresponding to a). The 
collapse parameter is a = 0.3. 



In addition to examining the temporal scaling of the 
pulse, we can also measure in detail the functional shape 
of the scaling solution. Here is it possible to distinguish 
between different subdiffusive processes. We have inves- 
tigated two candidate models for radial core spreading; 
the fractional diffusion equation (FDE) and the porous 
medium equation (PME). The fractional diffusion equa- 
tion is 
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where j = 2a denotes the order of a fractional time 
derivative 'io', ^2^. Solutions of this linear equation have 
the property that the width of a narrow pulse initial con- 
dition grows as t", where a < 1/2. If a = 1/2, the 
solution reduces to normal Fick diffusion. This FDE 
model is often used to describe processes which occur in 
spaces where there are temporal or spatial constraints, 
such as the flow of tracers through porous media E^ . 
The FDE has an analytic series solution in terms of Fox's 
H- Functions [2^ , which forms the self-similar scaling 
solution. We also examined the porous medium equation 
(PME), 



(2) 



This nonlinear model describes the spreading of a com- 
pact mound, and has the property that for a narrow pulse 
initial condition, the width grows as t^/^, and the scaling 
solution has a parabolic profile |2^. 

We fit radial core concentration data collapsed with a 
as a free parameter to the series solution of the FDE, 
and data collapsed with a = 1/3 to a numerical solution 
of the PME, as shown in figure 0^ and b respectively. 
We find that while both solutions model the collapsed 
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concentration profiles reasonably well within experimen- 
tal error, the PME has a smaller systematic discrepancy, 
since the profiles are better described as parabolic. The 
FDE solution has exponential wings and inflection points 
that are not obvious in the data. We note, however, that 
our projection visualization technique may simply be too 
insensitive to detect these tails. 

We also fit the non-segregating self-diffusion of the 
large grains to both models and find that the FDE gives 
a qualitatively better fit because in this case, the concen- 
tration profiles have tails within experimental resolution, 
while the parabolic PME solution does not. Examples 
of these fits to collapsed concentration profiles of mixing 
salt grains are shown in figures ^-d. In all cases, how- 
ever, fits to ordinary Pick diffusion with a = 1/2 are very 
poor. 

In conclusion, our results show that the axial transport 
of grains in a rotating tube is a subdiffusive process. This 
is true of both small particles comprising a segregated ra- 
dial core as well as for surface mixing of larger grains. In 
all cases, we find temporally self-similar concentration 
profiles that scale approximately as t^l'^ . These conclu- 
sions suggest that spontaneous axial segregation patterns 
in such tubes are more weakly damped, in the sense that 
they are sustained against slower mixing processes, than 
has been previously supposed. The goal of our future 
work is to elucidate the connection between axial band 
formation and the axial transport of grains, which is still 
unclear. 
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FIG. 4: (Color online) Collapsed concentration profiles of a 
black sand radial core in salt grains fit to: a) the fractional dif- 
fusion equation (yellow line), b) the porous medium equation 
(yellow line). Collapsed concentration profiles of mixing dyed 
black salt grains fit to: c) the fractional diffusion equation 
(yellow line), d) the porous medium equation (yellow line). 
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